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The condition for a lattice Dirac operator D to reproduce correct chiral anomaly at each site of a finite lattice for 
smooth background gauge fields is that D possesses exact zero modes satisfying the Atiyah-Singer index theorem. 
This is also the necessary condition for D to have correct fermion determinant ( ratio ) which plays the important 
role of incorporating dynamical fermions in the functional integral. 



1. INTRODUCTION 

Chiral anomaly plays an important role in 
particle physics. In QCD, the chiral anomaly 
breaks the global chiral symmetry of massless 
QCD through the internal quark loop. When the 
axial current couples to two external photons, the 
axial anomaly W provides a proper account of the 
decay rate of ir 77. For the flavor singlet ax- 
ial current coupling to two gluons, the presence of 
chiral anomaly but no associated Goldstone bo- 
son posed the problem which was resolved 
by 't Hooft, after taking into account of the topo- 
logically nontrivial gauge configurations, the in- 
stantons This may also provide an explana- 
tion why the rj { rj' ) particle is much heavier than 
the tt's ( K's ). 

Lattice QCD should be designed to provide 
nonperturbative and quantitive answers to all 
these and related problems. However, these 
goals could not be attained if the lattice Dirac 
fermion operator does not reproduce correct chi- 
ral anomaly or the fermion determinant ( ratio ) 
on a finite lattice without fine-tuning any param- 
eters. 

A basic requirement for lattice Dirac operator 
D is that, on a finite lattice with any prescribed 
smooth background gauge field which has integer 
topological charge, D possesses exact zero modes 
satisfying the Atiyah-Singer index theorem ( i.e., 
D is topologically proper ) without fine tuning any 
parameters. 



For any topologically proper D, the sum of the 
anomaly function over all sites must be correct, 
and is equal to two times of the topological charge 
of the background gauge field. If the anomaly 
function of D at some of the sites do not agree 
with the Chern-Pontryagin density, we can per- 
form the following topologically invariant trans- 
formation 101 



T{R): D' ^T{R)[D]= D{1 + RD)- 



(1) 
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with some operator R such that D' is local, then 
the anomaly function of D' would be in good 
agreement with the Chern-Pontryagin density at 
each site0. It suffices to choose R to be local in 
the position space and trivial in the Dirac space. 
Note that it is not necessary to perform any fine 
tunings of R since any local D' would give the 
correct chiral anomaly. The set of transforma- 
tions, {T(i?)}, form an abelian group with the 
group parameter space {i?} From the physi- 
cal point of view, we must require the existence of 
a chirally symmetric Dc — T{Rc)\D\ such that in 
the free fermion limit Dc{p) ijp,Pp, as p ^ 0, 
otherwise, D is irrelevant to massless QCD. In 
general, we assume that there exists 

Dc = D{I + RcD)'^ 

such that I?c75 + JaDc = 0. Then one obtains 

Dj5 + I5D - -2DRrj5D (2) 

This is the rejuvenated Ginsparg- Wilson rela- 
tion H]. In general, for any D, if the chirally 

^ Here we have assumed that the size of the finite lattice 
is large enough such that the finite size effects can be ne- 
glected, i.e., the size of the lattice is much larger than the 
localization length of D' . 
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symmetric Dc = T{Rc)[D\ exists, then D must 
satisfy the Ginsparg- Wilson relation. Thus the 
Ginsparg- Wilson relation does not specify the 
topological characteristics of a Dirac operator. 
Only a topologically proper D ( or its transform 
D' = D{\ + RD)-^ ) is suitable for lattice QCD. 

2. THE ANOMALY FUNCTION 

Consider a lattice Dirac operator D which 
breaks the usual chiral symmetry according to 
1^75 + 75-D = B, where i? is a generic irrelevant 
operator, and the lattice is finite with periodic 
boundary conditions. Then we can obtain Q 



which satisfies the conservation law, 



(3) 



If D does not possess exact zero modes in the 
background gauge field, then the fermionic aver- 
age of d^J^{x) is 

{d^Jl{x)) - i tr [{BD-^){x,x) + {D-'B){x,x)\ 



= An{x) : anomaly function 

However, for any D which does not possess exact 
zero modes in the background gauge field, the 
sum of the anomaly function over all sites must 
vanish due to the conservation law, 

J2{d'^4{x))=J2Ar.ix)^0 



This implies that if An (x) is not zero identically 
for all X, then it must fluctuate from positive to 
negative values with respect to x. The latter case 
is exactly what happens to the anomaly function 
of the standard Wilson-Dirac fermion operator. 

On the other hand, if D possesses exact 
zero modes in topologically nontrivial background 
gauge fields, then is not well defined. In this 
case, one introduces an infinitesimal mass m and 
replace D hy D = D + mf[D] ( where f[D] is 



any functional of D, which has eigenvalue one for 
the exact zero modes oi D ), and finally take the 
limit ( m ^ ), then we obtain 



A'4- 



(d'^^ix)) 



An{x)+2j2[cl>t{x)]Ut{^) 

5-1 

-2^[</.r(x)]t0-(a;) (4) 



where 0+ and (j)^ are normalized eigenfunctions 
of D with eigenvalues As = At = and chiral- 
ities +1 and —1 respectively^, and the anomaly 
function is 



An{x) = lim i tr 



{Bb-^){x,x) + {t)-^B){x,x) 



3. THE INDEX THEOREM 

Summing Eq. (Q) over all sites and using the 
conservation law Eq. (^, we obtain 



(5) 



This is the index theorem for lattice Dirac oper- 
ator on a finite lattice, in any background gauge 
field. However, it does not necessarily imply the 
existence of exact zero modes nor the compliance 
with the Atiyah-Singer index theorem. 

Since D does not have exact zero modes in a 
trivial gauge background, the index of D must 
be proportional to Q of the smooth background 
gauge field. If Q is an integer, then the pro- 
portional constant must be an integer, otherwise 
their product in general cannot be an integer. De- 
noting this integer multiplier by c[D\, we have 



N^- c[D] Q uf 



(6) 



where Uf is number of fermion flavors. Here we 
have assumed that c[D] is constant for smooth 
background gauge fields. This is a reasonable 
assumption since c[D] is an intrinsic character- 
istics of D. However, when the gauge field be- 
comes rough, we expect that the index theorem 
with integer c[D] would break down. If one insists 



^ Here we assume that D is normal, DD^ = D^D, and D 
satisfies the hermiticity condition = 75D75. 
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that it holds even for rough gauge configurations, 
then c[D] cannot be an integer constant due to 
the highly nonlinear effects of the gauge field. In 
general, c[D] is a rational number functional of 
D, which in turn depends on the gauge configu- 
ration, but it becomes an integer constant only 
for smooth gauge configurations. The topological 
characteristics of D, c[D], was first discussed in 
ref. ^ and was investigated in ref. for the 
Neuberger operator l^. 

4. A GENERAL SOLUTION FOR An{x) 

Now we try to obtain a general solution of 
An{x) satisfying (||) and (|6|), i.e., 

i ^N_-N+= c[D] Q Uf. 

X 

Consider the gauge configuration with constant 
field tensors, 

L1L2 J^3^4 

and other F^s are zero, where qi and q2 are inte- 
gers. Then the topological charge of this config- 
uration is 

X 

= n qi q2 

which is an integer. If D is local, then An{x) is 
constant for all x. We immediately obtain 

A-n {x^ 

= ^c[D]e^,x^F;,{x)Fl{x + A + i>)tr{t"t''} 

This is the particular solution of An (x) . The com- 
plementary solution to the homogeneous equation 
X^a; -^n{x) = in general can be written as 

^("Hx) = [^A'(^) -G^x- A)] = d^G^ix) 



where G^{x) is any local and gauge invari- 
ant function. Then the general solution of the 
anomaly function is the sum of the complemen- 
tary solution and the particular solution. 



An{x) 

= ^c[D]e^,^,F^,{x)Fl{x + A + i^Mt^t"} 
+ Y[G^ix)-G^{x-f,)] (7) 

In general, G^{x) depends on the background 
gauge field and D. Only when 9pG^t(x) vanishes 
identically for any smooth gauge background, D 
can reproduce the correct chiral anomaly pro- 
vided c[D] = 1. ( e.g., the Neuberger operator 
1^ ). On the other hand, the pathologies of the 
standard Wilson-Dirac operator are c[D] — and 
A''^\x)^Q. 

Now we consider D with d^G^{x) = in a 
constant gauge background. Then we introduce 
local fluctuations to the background gauge field 
and/or increase its topological charge Q. Then 
c[D\ may remain constant provided that the lo- 
cal fluctuations of the gauge fields are not too vi- 
olent and/or Q is not too large. When the gauge 
background becomes so rough that D is non-local, 
then d^G^{x) ^ but the index theorem may still 
hold with the same c[D\. If we keep on increasing 
the roughness of the gauge background, then D 
will undergo a topological phase transition, and 
c[D\ will become another integer or even a frac- 
tion 0. 
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